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Abstract. The limited number of active pixels in a Geiger-mode Avalanche 
Photodiode (G-APD) results not only in a non-linearity but also in an additional 
fluctuation of its response. Both these effects are taken into account to calculate the 
amplitude resolution of an ideal G-APD, which is shown to be finite. As one of the 
consequences, the energy resolution of a scintillation detector based on a G-APD is 
shown to be limited to some minimum value defined by the number of pixels in the 
G-APD. 
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Multipixel Geiger-mode avalanche photodiodes (G-APDs) pQ are solid state 
photodetectors with a fast growing field of applications. A G-APD can be thought of as 
a matrix of independent photon counters (pixels / cells) connected to a common load. 
A photoelectron produced in the G-APD volume by a photon can trigger an avalanche 
breakdown in one of the cells (fire the cell). To an input signal of n photoelectrons 
the G-APD responds with an output signal of N fired cells, the dependence N = f(n) 
saturates due to the limited number m of the cells in the device, and the fact that more 
than one photon can enter one cell without changing the response of that cell. The 
non-linearity of a G-APD response is an established experimental fact (see, e.g [2]). A 
correlation between saturation of the G-APD signal and an increase of its dispersion 
was also observed [2]. 

In this work we calculate the amplitude resolution of an ideal G-APD (which we call 
the "intrinsic" resolution) and show that it is limited (non zero) due to an additional 
statistical noise in the process of distributing n photoelectrons over m cells. By an 
ideal G-APD we understand a device with the excess noise factor [3] equal to unity. In 
particular, we assume that there is no crosstalk between cells (so it is never the case 
that one photon triggers more than one cell) and that each photon incident hits some 
cell. 

The considered problem is equivalent to a well-known problem in mathematical 
statistics of distributing (randomly) n balls (photoelectrons) into m urns (cells), see 
e.g. [I]. The number N of urns containing one or more balls is a random variable, its 
expected value and variance are: 

N = m [1 - (1 - mr x ) n \ 

(1) 

a% = m (m — 1) (1 — 2 m 1 ) n + m (1 — m 1 ) n — m 2 (1 — m 1 ) 2n 

The distribution of N is approximately normal when m, n — > oo and the ratio a = n/m 
is bounded pE]: 

N = m(l-e~ a ) 

(2) 

a 2 N = me- a [l - (l + a)e" Q ] . 

In practice the number of cells in a G-APD is usually greater than ~ 100, which 
justifies using the asymptotic formulae ([2]) in the following analysis. 

The amplitude resolution of a detector is a measure of its ability to distinguish 
between input signals having different amplitudes j5] . For a device with a linear response 
the resolution can be shown to be: R = 2.35 ■ o~/A, where A is the amplitude of the 
detector response, and A and a are the mean value and standard deviation characterizing 
the distribution of A (which is assumed to be normal). In case of a non- linear response, 
the formula has to be modified. The scheme of our consideration is the following: 

A 1 M A:(A,a) ^ A 2 :(A 2 ,a 2 ). 
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To an input signal Ai (constant value) the detector responds with a signal of amplitude 
A. Due to fluctuations induced by the detector, A is statistically distributed. If the 
function / relating A to A\ is known, one can reconstruct an estimate of the input signal 
from the mean detector response. The restored amplitude A 2 is no longer a constant, 
like A\, but is also a random value, with mean and standard deviations (at a< A): 

A 3 = f- 1 (A)=A 1 , a 2 =(f- 1 )' A -a. (3) 

The parameters of the restored distribution, rather than the measured one, define 
the amplitude resolution of a detector: 

R = 2.35 ■ (a 2 / A 2 ) . (4) 

The theoretical function N — f(n) is known (j2J), and this allows us, using relations 
(IHJ), to reconstruct the signal at the G-APD input (A 2 is the reconstructed amplitude) 
from the measured distribution of iV (in ([3]) A = N, a = a^): 



a 



A 2 = n , a 2 = a N e 

By substituting A 2 and a 2 in (j3J) for the intrinsic amplitude resolution R of a 
G-APD we get: 

Ro = 2.35 • m- 1 / 2 ■ $ (a) , 

(5) 

$ ( a ) = a- 1 (e a - 1 - a) 1 / 2 . 

The function <&o(oe) is plotted in Figure 1. For any given value of m the value of 
Rq increases (the resolution becomes worse) with the ratio n/m = a increasing. 




Figure 1. The functions $0(0;) and <&(a) defined in j5|) and ([6|). 



When a G-APD is used in a scintillation detector the intrinsic resolution of the G- 
APD should be summed quadratically (this we prove in Monte-Carlo simulations below) 
with the resolution determined by the statistics of photoelectrons -R s t a t = (the 
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mean number of photoelectrons n relates to the number of photons as n = PDF ■ n p h, 
where PDE is the photon detection efficiency |3]). Note, that here we assume that the 
energy resolution of the scintillator [6] is equal to zero. The energy resolution of such 
an ideal scintillation detector is then: 

R = 2.35 -m- 1 ' 2 -^(o) , 

(6) 

$(«) = (u- l + <$>l(u)) 1/2 . 

The function $(a) (see Figure 1) has minimum value of $ m i n = 1.2426 at a = 
1.5936. Accordingly, the minimum value for the energy resolution of an ideal G-APD 
based scintillation detector (taking into account only the statistics of photoelectrons) is 
R min = 2.92/ y/m. 

The obtained result on the energy resolution was verified in Monte-Carlo 
simulations. In the simulations the energy resolution of an ideal scintillation detector 
was calculated as the follows: a) a number of photoelectrons n was randomly generated 
according to Poisson statistics; b) those n photoelectrons were randomly (uniformly) 
distributed over m cells resulting in N occupied (fired) cells, c) from the distribution 
of N the number of photoelectrons (the energy) and its standard deviation were 
reconstructed using ([3]), and the energy resolution calculated according to (jlj). As is 
seen in Figure 2, the energy resolution obtained in simulations is in agreement with the 
values calculated according with (J6]). 

Physically, while a larger mean number of photons implies a smaller fluctuation 
in that mean, at some point saturation effects dominate (multiple photons enter single 
cells) and resolution degrades. This means that some care should be taken in each 
application to optimize design taking this into account. 

Summary 

The variation of a G-APD response induced by the statistics of distributing n 
photoelectrons over m G-APD cells causes its amplitude resolution to degrade as the 
ratio n/m increases. This implies, for example, in that the energy resolution of a 
scintillation detector using a G-APD is limited to -R m i n = 2.92 / ^fm. 
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Figure 2. (Top) Simulated values (dots) of the energy resolution of an ideal G- 
APD based scintillation detector as compared with the values (solid line) calculated 
according with ([6]). The dashed lines show two contributions to the energy resolution: 
the statistical term (l/^/ma) and the intrinsic resolution of the G-APD (Rq). 
(Bottom) the same for G-APDs with different number of cells m. 



